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We investigate the effect of the experimentally observed Jahn- Teller distortion of the oxygen oc- 
tahedra in LaTiOs on the magnetic exchange. We present a localized model for the effective hopping 
between nearest-neighbor Ti ions and the intra-site Coulomb interactions, based on a non-degenerate 
orbital ground state due to the static crystal field. The latter corresponds to an orbital order which 
recently has been confirmed experimentally. Using perturbation theory we calculate, in addition to 
the Heisenberg coupling, antisymmetric and symmetric anisotropy terms of the superexchange spin 
Hamiltonian which are caused by the spin-orbit interaction. Employing our superexchange Hamil- 
tonian, we deduce the magnetic order at low temperatures which is found to be in good agreement 
with experiment. 

PACS numbers: 71.10.-w, 71.27.+a, 71.70.Ch, 75.10.Dg, 75.25.-fz, 75.30.Et 



I. INTRODUCTION 

In the seventies, the orthorhombic perovskite LaTi03 
was considered as a typical example of an antiferromag- 
netic Mott insulator (Tn = 146 K). 1 The ground state 
of the Ti ion is trivalent with a single electron in the 
d shell. Recently, LaTiOs has attracted attention when 
an ordered magnetic moment of 0.46 fiB was reported. 2 
This value is surprisingly small for a single electron with 
quenched orbital moment, for which one would have ex- 
pected 1//B- In LaTiC>3, this value may be reduced 
by about 17% due to quantum fluctuations of the 3D 
Heisenberg model, and by about further 14 % due to the 
on-site spin-orbit coupling, (in conjunction with the crys- 
tal field, see below), leading to an overall estimate of 
0.72 fj, B . 

A previous attempt to explain this unexpected finding 
has neglected the Jahn- Teller distortion of the oxygen 
octahedra in LaTi03 and assumed that the symmetry of 
the unit cell is strictly cubic. In such a situation the ti g 
ground state of the Ti ion is three-fold degenerate and 
the orbital moment is unquenched. Hence, it was pro- 
posed to consider LaTiC>3 as an orbital liquid in order to 
explain the reduction of the ordered moment by orbital 
fluctuations. 3 

However, recent experiments give a strong indication of 
the importance of the Jahn- Teller distortion in LaTiC>3, 
and in particular enable, using recent structural data, to 
estimate the splitting it induces in the ty, g levels: There 
is a crystal-field gap of about 0.2 eV between the non- 
degenerate ground state and the next excited level. This 
value has been confirmed by a study of photoelectron 
spectroscopy 5 and is at least one order of magnitude 
higher than any superexchange energy in LaTiOs. 6 Con- 
sequently, the orbital order at low temperatures is not in- 
duced by the superexchange. Rather, the orbital degree 
of freedom is frozen by the crystal field. The scenario of 



suppressed orbital fluctuations has also been confirmed 
by a recent LDA+DMFT study. 7 The assumption of Ref. 
8 that the crystal-field splitting, the superexchange, and 
the spin-orbit coupling are all of the same order is incon- 
sistent with the photoelectron spectroscopy of Ref. 5. In 
addition, a higher value, 0.57 /j,b, of the ordered moment 
has been recently reported, 4 making the discrepancy be- 
tween experiment and theory (which gives 0.72 /ib) even 
smaller. 

The Jahn- Teller effect in LaTi03 is caused by the 
twisting of the Ti-0 bonds with respect to each other 
(i.e. by differences between the 0-0 bond lengths), 
rather than by differences between the Ti-0 bond 
lengths. The non-degenerate ground-state orbital due 
to crystal-field calculations given in Ref. 4 is consistent 
with the orbital order found in NMR measurements of 
the Ti-3d quadrupole moment. 9 The presence of orbital 
order at low temperatures has been also concluded from 
measurements of the dielectric properties and the dy- 
namical conductivity. 10 As opposed to these findings, an 
orbital contribution to the specific heat, which is pre- 
dicted by the orbital-liquid model, has not been found in 
experiment. 11 

Hence, from the recent experiments it must be con- 
cluded that the orbital-liquid model is inappropriate for 
LaTiOs. Moreover, it has been proven, by exact sym- 
metry arguments, that due to a hidden symmetry the 
superexchange Hamiltonian used in Ref. 3 cannot repro- 
duce the observed magnetic order of LaTiOs. 12,13 

Another model to explain the magnetic properties of 
LaTiOs proposed the lifting of the ti g degeneracy by the 
crystal field resulting from the eight La ions which sur- 
round each TiOg octahedron — assuming undistorted oc- 
tahedra, i.e. neglecting the Jahn- Teller effect. 14 This 
model predicts a realistic non-degenerate ground-state 
orbital for each Ti ion, 6 ' 9 and yields plausible values for 
the Heisenberg couplings between nearest-neighbor Ti 
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ions. 



II. THE MODEL 



However, there are two points missing in the calcula- 
tions of Ref. 14. (i) The crystal field due to the eight 
nearest La ions, which gives a nearly equidistant split- 
ting scheme between the three ti g orbitals, is only a first 
approximation of the full crystal field due to all ions 
of the solid. It is preferable to treat the electrostatic 
crystal field more accurately, employing the Madclung 
sum. Such a treatment shows that the Jahn- Teller effect 
leads to a non-degenerate ti g ground-state orbital and 
two quasi-degenerate excited orbitals. 4 (ii) Terms of the 
exchange Hamiltonian which break the spin-rotational in- 
variance and cause the magnetic order have not been con- 
sidered. Hence, the origin of the observed magnetic order 
is not fully understood so far. 

In the present paper, we investigate a model for the 
magnetism of LaTiC>3, which starts from a point-charge 
calculation of the static crystal field for the Ti ions via a 
full Madelung sum over the crystal (as was already dis- 
cussed in Ref. 4). Taking into account the recent struc- 
tural low-temperature data and using a Slater-Koster 
parametrization of the Ti-0 hopping, we calculate an ef- 
fective hopping matrix between the d orbitals of nearest- 
neighbor Ti ions. Treating this Ti-Ti hopping and the 
on-site spin-orbit coupling as perturbations, we calculate 
the superexchange coupling between the non-degenerate 
crystal-field ground states of the Ti 3+ ions. In treat- 
ing the Ti 2+ ions, which appear as intermediate states 
of the exchange processes, we take into account the full 
on-site Coulomb correlations in terms of Slater integrals 
and diagonalize the Coulomb Hamiltonian together with 
the crystal-field one. The spin-orbit coupling gives rise 
to antisymmetric and symmetric anisotropies of the spin 
Hamiltonian. We calculate the isotropic part of the ex- 
change coupling and both kinds of the anisotropies to 
leading orders. Using our exchange Hamiltonian, we de- 
termine the classical ground state which gives the direc- 
tions of the spins in the ordered phase. The experimental 
data reveals a G-type antiferromagnetic order along the 
crystallographic a axis, which is accompanied by a small 
ferromagnetic moment along the c axis. 4 Our calculation 
reproduces this order. In addition, we find a small A- 
type moment along the b axis, which has not yet been 
detected experimentally. 

In the next section we present the details of our model. 
Section III is devoted to the perturbation expansion 
yielding the microscopic spin Hamiltonian. Section IV 
discusses the macroscopic magnetic Hamiltonian and the 
resulting magnetic order of the classical ground state. It 
includes as well a detailed comparison with existing ex- 
perimental data. Finally, we summarize our results in 
Sec. V. 



A. The crystal field 

The unit cell of LaTiC>3 contains four Ti ions, see Fig. 
1 and Table I (the Tables are given on pages 18ff.), and 
has the symmetry of the space group Pbnm (No. 62 in 
Ref. 15). The symmetries of this space group are listed 
in Table II. Given the position of one La, Ti, 01, and 02 
ion each (see Table I), the positions of all other ions in 
the unit cell follow from the space-group symmetries. In 
order to use these symmetries conveniently, we employ 
in our calculation the orthorhombic orthonormal 16 basis 
for the Ti-d orbitals 

\xy), |2z 2 ), \yz), \xz), \x 2 - y 2 ), (1) 




FIG. 1. The crystallographic structure of LaTiOs. The 
ten Ti ions, which constitute the twelve inequivalent near- 
est-neighbor Ti-Ti bonds are enumerated. For simplicity, 
oxygen octahedra are only shown around the four crystal- 
lographically inequivalent Ti sites. La ions from two layers 
are shown as small spheres. For example, the sites 2 and 6 
are crystallographically equivalent but the bond 61 emerges 
from the 12 bond by a glide reflection, so that the effective 
hopping matrix for the bond 16 is different than that of the 
12-bond: It is the transposed one, see Table IV. 

where the x,y and z axes correspond to the crystallo- 
graphic a, b and c axes. In a cubic perovskite, the first 
two orbitals would correspond to the e g orbitals and the 
three others to the t 2g orbitals. [Note that the pseudo- 
cubic basis for the d orbitals, which is frequently used, is 
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obtained from Eq. (1) upon rotating the x and y axes by 
45° around the z direction.] 

Using the structural data of Ref. 4 (taken at T = 8 K), 
we have calculated the spectrum and the eigenstates 
for the Ti ion located at (0,1/2,0), employing a point- 
charge calculation of the static crystal-field Hamiltonian. 
This calculation uses the full Madelung sum over the 
crystal (which is evaluated as an Ewald sum, see Ap- 
pendix A). It requires the second moment, (r 2 ), and 
the fourth moment, (r 4 ), of the effective ionic radius of 
the Ti 3+ -ion, see Appendix A. We have used the values 

(r 2 ) = 0.530 A 2 and (r 4 ) = 0.554 A 4 . 17 The results of 
the crystal-field calculation, which are listed in Table III, 
exhibit a typical Jahn- Teller ti g splitting scheme, where 
a non-degenerate ground state is clearly separated from 
the two quasi-degenerate excited states. 

The orbital order due to the static crystal field is shown 
in Fig. 2. The ground-state orbital is given by 



|0) = 0.770|yz) ± 0.636|x 2 - y 2 ). 



(2) 




This state has approximately the 2z 2 structure in the 
coordinate system in which the y and z axes are ro- 
tated by ±56° around the x axis. 4 The relative sign of 
the linear combination alternates between neighboring ab 
planes according to the mirror planes at z — 1/4, etc. In 
the pseudo-spin language, 18 we have ferro-orbital order in 
the ab planes and canted antifcrro-orbital order between 
the planes. We note that this ground state is in perfect 
agreement with experiment. 9 The ground state cited in 
Refs. 9 and 14 is given, to a good approximation, by 




0.816 0.577 



It practically coincides with our ground state, 

|(0|0')| 2 = 99.06%. (4) 



The ground state |0) [see Eq. (2)] of the crystal field, 
which is occupied at each Ti site by a single electron, is 
the starting point of our model. The perturbative cal- 
culation outlined below is employed in order to evaluate 
the magnetic superexchange coupling between Ti ions in 
this state. 



FIG. 2. The orbital order of the Ti ions resulting from the 
calculated crystal field (the energy scale of the superexchange 
which could in principle also affect the orbital degree of free- 
dom is about one order of magnitude smaller than the crys- 
tal-field gap). The order is ferro-orbital in the ab planes and 
canted antiferro-orbital between the planes. 



B. The Hamiltonian 

The perturbation scheme is applied to a two-site clus- 
ter consisting of two nearest-neighbor Ti ions, denoted 
by m and n. The unperturbed Hamiltonian acting on 
such a cluster is given by 

where H^ in is the static crystal field and H£>„ describes 
the intra-ionic Coulomb correlations of a doubly occupied 
d shell. The perturbation calculation requires a selected 
set of eigenstates and the corresponding eigenenergies of 
H® nn . In our case, the eigenstates span a Hilbert space 
which consists of two sectors. In the first, termed the 
Ti 3 + sector, both Ti ions are trivalent. In the second, 
called the Ti 2+ sector, one of the Ti ions is divalent (two 
d electrons on the same site) and the other is four-valent 
(an empty d shell). In the ground state of H^ n , which 
belongs to the Ti 3+ sector, both Ti ions are in the one- 
particle ground state of H^ n , modulo spin up or down 
on each site. This leads to a four-fold degeneracy of the 
ground state of the cluster. The complete Ti 3+ sector 
has 100 basis states (made up of five orbital and two 
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spin states on each of the two sites) . The Ti 2+ sector has 
90 basis states (30 spin-triplet states and 15 spin-singlet 
states of the doubly occupied d shell and a factor of 2 
because each of the Ti ions can be doubly occupied while 
the other has an empty d shell). 

In order to calculate the spectrum of H mn , we ap- 
ply H^ n on the Ti 3+ sector, and both H^ n and H^ n 
on the Ti 2+ sector. H^ in is parametrized in terms of 
the Slater integrals F% and -F4, 19 and the effective Ti- 
Ti charge-transfer energy U e g. This energy is the dif- 
ference between the four-fold degenerate ground state 
of the cluster (which is the lowest level of the Ti 3+ 
sector) and the lowest level of the Ti 2+ sector (where 
H^ n and H'^ nn are diagonalized simultaneously). Wc 
use F 2 = 8F4/5 = 8.3 eV from an atomic Hartree-Fock 
calculation, 20 and U c s — 3.5 eV from the analysis of the 
photoemission spectra and first-principles calculations. 21 

The perturbation part of the Hamiltonian, V mn , con- 
sists of an effective Ti-Ti tunnelling term, H^, and the 
on-site spin-orbit interaction, H^ n , 

v mn mn ' -^mn - v u / 

The tunnelling Hamiltonian is given in terms of an effec- 
tive hopping matrix, t mn , between the m and the n Ti 
ions, 

H^ = J2€J mi d nj +h.c, (7) 

ij 

where the indices i and j enumerate the eigen-orbitals of 
H^ n (see Table III). The spin-orbit coupling is given by 

H r s n °„ = A^l fe s fe , (8) 

k—m.n 

where lfc denotes the angular momentum operator of the 
Ti ion at the k site, its spin operator, and A the spin- 
orbit coupling strength. Wc use A = 18meV. 14 

The 14% reduction of the magnetic moment, alluded 
to in Sec. I above, is obtained upon diagonalizing to- 
gether H^[ n and H*° n for a single Ti 3+ ion. The reduction 
in the magnetic moment due to the not fully quenched 
orbital moment is estimated as follows. The spin along 
a certain selected direction does not commute with Hf° n 
and therefore it is not a good quantum number for the 
combined Hamiltonian ff^ f „ + -ff"„. However, the eigen- 
states of this combined Hamiltonian are symmetric or 
antisymmetric with respect to time-reversal. This leads 
to five Kramers doublets for the single Ti 3+ ion. Since 
the ordered magnetic moment is mainly of the G-type, 
and is oriented along the x axis, 4 we use those doublets 
to find the expectation values of the angular momen- 
tum. By choosing the largest possible polarization of the 
magnetic moment along the x axis out of all the linear 
combinations of the ground-state doublet, we find the ex- 
pectation value (1% + 2s£) /i B = 0.86 /x B . This effect is 
not included in our perturbation calculation of the mag- 
netic exchange. However, it does explain partially why 



the observed ordered moment along the x axis is reduced 
with respect to 1 u^. 

The dominant hopping process between two nearest- 
neighbor Ti ions is mediated via the oxygen ion which is 
nearest to both of them. Let t l £ be the hopping matrix 
element of an electron in the p orbital a on the oxygen ion 
into the i state of the Ti ion located at m. The effective 
hopping between the Ti ions is then given by 

Here, A e g- is the charge-transfer energy, which is required 
to put an electron from an O ion on a Ti ion, and a de- 
notes one of the three p orbitals on the oxygen (in or- 
thorhombic coordinates), 

|z>,|y>,|*>. (10) 

Using the structural data taken at T = 8K, 4 to- 
gether with elementary geometric considerations, the Ti- 
O hopping matrix elements can be expressed in terms 
of the Slater-Koster parameters V p da and V p dn- 22 Wc 
use the values V p da = — 2.4eV, V v dn = 1.3 eV, and 
A ff = 5.5 eV, 14 ' 21 in conjunction with Eq. (9) to com- 
pute the effective hopping matrices pertaining to the unit 
cell. The results are listed in Table IV, which also gives 
the symmetry properties of the hopping matrices between 
different Ti-Ti bonds. The four incquivalent Ti sites of 
the unit cell form twelve nearest-neighbor Ti-Ti bonds 
which are crystallographically inequivalent, i. c. they do 
not evolve from each other by Bravais translations. These 
bonds connect the ten Ti ions indicated in Fig. 1. By the 
symmetry operations of the space group Pbnm, the eight 
effective hopping matrices between Ti ions belonging to 
the same ab plane and the four matrices for Ti-Ti bonds 
along the c direction, respectively, can be expressed by a 
single matrix each. For example, all twelve hopping ma- 
trices are given by the two matrices for the Ti-Ti bonds 
mn = 12 (planar) and mn — 13 (inter-planar), respec- 
tively. 

Despite of the different orbital ordering within the ab 
planes and between them, the hopping amplitudes be- 
tween the crystal-field ground states in and between the 
planes are of the same order of magnitude, i. e. roughly 
1*12 1 ~ |*i3 1- In strictly cubic symmetry, those are equal 
(see Appendix B). The rotation of the oxygen octahedra 
around the c axis, the tilting, and the distortion cause 
some difference between \tj^\ and | *?3 1 - 

For convenience, we present in Table V an overview of 
the parameters used in our calculation. 



C. The Ti-O hybridization 

Our model does not include the covalent contribution 
to the crystal field, arising from hybridization between 
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the Ti-3ci and 0-2p states. This mechanism mixes ex- 
cited states of the static crystal-field into the Ti 3+ ground 
state, i.e., there is an admixture of Ti 2+ states accom- 
panied by an admixture of holes on the oxygen sites. 

Following Ref. 21, we now estimate the effect of the 
pd hybridization. When that hybridization is absent, the 
effective parameter U e g defines the energy difference be- 
tween the ground state of the Ti 3+ sector and the lowest 
state of the Ti 2+ sector in a two-site cluster consisting 
of two Ti ions. When the pd hybridization is present, 
these two types of d states correspond to two bands, 
from which two pd hybridized bands evolve according 
to the covalent crystal field. These hybridized bands 
have, in general, considerable dispersion: Their peak- 
to-peak separation, which is seen in the combined pho- 
tocmission and inverse photocmission spectra, is given 
by the band gap £ ; gap =1.6eV, and the distance between 
the band edges is given by the optical gap -E op t=0.2eV, 
which is experimentally observed as the Mott gap. The 
mean bandwidth between the two pd hybridized bands is 
then W = E eap - E opt = 1.4 eV. 

Since the bands are rather dispersive, the question 
arises whether a localized picture is suitable to describe, 
even approximately, the LaTiC>3 system. In order to 
study this point, we have analyzed the covalent crystal 
field of a cluster consisting of a single Ti ion, and the six 
oxygen ions predominantly hybridized with it (the calcu- 
lation has been carried out for Ti No. 1 in Fig. 1). This 
is accomplished by diagonalizing the Hamiltonian 

H pd = H cf + H c + H£T, (11) 

for a Ti0 6 -cluster. Here H ct describes the static crys- 
tal field, H c is the Coulomb interaction, and H^ n is the 
pd-tunnelling, 

^ n = E^4p„ tt +h.c. (12) 

nioi 

In Eq. (12), p na destroys an electron on the n-th oxygen 
site in the a-orbital, given in Eq. (10). As in the calcula- 
tion of the Ti-Ti hopping amplitudes, the pd hopping am- 
plitudes, t l {^ , are expressed in terms of the Slater-Koster 
parameters V p da and Vpdn, using the structural data of 
Ref. 4. The entire space of the basis states of the TiO@- 
cluster consists of a Ti 3+ sector where the p orbitals are 
all occupied, and a Ti 2+ sector where there is a hole in 
one of the p orbitals. The cigenstates of the Hamiltonian 
(11) have the form 

|V>) = V2~ ^d 1 ) + V^d~l\d 2 ), (13) 

where rid is the occupation number of the Ti-cZ shell 
(1 < rid < 2), jd 1 ) is a state with a single electron in 
the d shell and fully occupied p shells on the surround- 
ing oxygen ions, and \d 2 ^ is a state with two electrons in 
the d shell and a hole in the p shell of one of the oxygen 
ions. We find that in the ground state rid — 1-334 , i. e., 



there is a p hole on one of the neighboring oxygens with 
probability of 33.4%. 

This calculation allows for the analysis of the eigen- 
states of the combined static and covalent crystal fields. 
Projecting the five lowest eigenstates onto the Ti 3+ sec- 
tor, (which corresponds to the states jd 1 )), gives to a 
very good approximation the same eigenstates as for the 
static crystal field alone, see Table VI. This finding ex- 
plains why, despite the admixture of Ti 2+ states |d 2 ), 
the agreement with the experiment remains perfect, as 
evidenced in Eq. (4). Indeed, the experiment measures 
the Ti 3+ part, jd 1 ), of the combined static and covalent 
crystal field, and apparently is not sensitive to the Ti 2+ 
admixture \d 2 ). Table VI also shows that the ti g split- 
ting remains almost the same as in the absence of the 
covalent contribution, whereas the distance between the 
tig and e g energies is enhanced. 

Since it is extremely complicated to include in the mag- 
netic superexchange calculation the hopping between the 
pd hybridized states, we choose to consider the hopping 
between the Ti 3+ states only. The results listed in Table 
VI, which show that the projections of the eigenstates 
of the combined static and covalent crystal fields onto 
the Ti 3+ sector are almost the same as in the static-only 
case, ensure that the Ti 3+ states we use are an appropri- 
ate starting point for the superexchange calculation. 



D. Comparison with other models 

It is instructive at this point to dwell on several dif- 
ferences between our model and three other calculations, 
reported in Refs. 7, 14, and 23. 

(i) As is mentioned above, Ref. 14 takes into account 
the crystal field due to the La ions only. In addition, 
the intra-ionic Coulomb correlations are approximated 
according to the Kanamori scheme, which ignores the 
splitting of the spin-triplet states. As opposed to this cal- 
culation, we take into account the crystal field for both 
the Ti 3+ and the Ti 2+ configurations, and employ the full 
intra-ionic Coulomb correlations for the latter. The last 
point is particularly important: The spin-triplet states 
as intermediate Ti 2+ states cause a ferromagnetic cou- 
pling while the spin-singlet states induce an antifcrro- 
magnetic coupling, leading to a competition between fer- 
romagnetic and antiferromagnetic contributions to the 
magnetic exchange. On the other hand, we omit the 
small pp hybridization between the oxygen states, which 
was included in the calculation of Ref. 14. 

(ii) The LDA+DMFT calculation 7 gives a ground 
state (denoted here |0")) whose projection on the 
experimentally-deduced ground state, Eq. (3), is 

|(0'|0")| 2 = 87.8%, whereas we find 99.06%, see Eq. 
(4). An even larger difference, (which is partially ex- 
plained by the difference in the ground-state orbitals) is 
found between the nearest-neighbor hopping amplitudes, 
coupling the Ti ions in that ground state: The values 



5 



cited in Ref. 7 are about half of the ones we use, with 
the in-plane amplitude being slightly smaller than the 
inter-plane one. 

(iii) A recent calculation of the crystal field at room 
temperature, including the covalency contribution and 
the spin-orbit coupling, 23 has yielded basically the same 
t2g splitting scheme as ours (due to Table VI), whereas 
the spacing between the ti g and e g states turned out to 
be bigger than in our calculation, about 0.9 eV. Analo- 
gously to the way we determined the reduction of the 
magnetic moment due to the spin-orbit coupling in Sec. 
II B, the ground state found in Ref. 23 has entangled 
spin-up and spin-down states, i. e. the orbital part is not 
separable from the spin part. The ground state is given 
there in the fashion which has the largest possible mag- 
netic polarization along the quantization axis. Denoting 
the spin- up part of this ground state by |0'"), it turns 
out that the squared overlap with the experimentally de- 
termined orbital is |(0'|0"') | 2 = 92.47%. In Ref. 23 the 
reduction of the G-type moment due to the crystal field 
and the spin-orbit coupling is found to be 9.5 %, whereas 
we find 14%, see Sec. II B. 

We will continue the comparison with other models 
when we discuss results of our calculation in Sec. IV. 



III. THE EFFECTIVE SPIN HAMILTONIAN 

Our aim is to obtain from the full Hamiltonian, H mn — 
Hmn + Vmn' an effective spin Hamiltonian, h mn , which 
acts within the Hilbert space of the four-fold degenerate 
ground state of the unperturbed Hamiltonian H mn . 

In general, an operator which acts in the ground-state 
space of the two Ti ions located at sites m and n, consists 
of linear combinations of following terms, 



d\nQa' ^nOcr' d n 0a n ^mOcr > 



(14) 



where d]j 0(7 (d n0r7n ) creates (destroys) an electron in the 
crystal-field ground state at site n, of spin component a n . 
Since there is a single electron at each Ti site, the cre- 
ation and annihilation operators can be written in terms 
of site spin- 1/2 operators, S„, 



1_ QZ 

2 D n- 



(15) 



Any operator acting within the ground-state space of 
the two Ti ions can be represented in terms of the 16 
operators 



1 (constant), 
S% (single-ion anisotropics), 
S m Sf^ (inter-site spin couplings) 



(16) 



anisotropies, and consequently the effective spin Hamil- 
tonian (omitting constant terms) takes the form 

h>mn — S m ^4 mn S n , (f^) 

where A mn y = A* lm ) is the 3x3 supcrcxchangc ma- 
trix. This matrix may be decomposed into a symmetric 
part and an antisymmetric one. The three components of 
the latter constitute the Moriya vector D mn (= — D nm ). 
Extracting further the isotropic part of A mn , i.e., the 
Heisenberg coupling J mn , the effective spin Hamiltonian 
is cast into the form 

hmn = Jmn^m^n + D m „(S m X S„) + S rn A mn S n . (18) 

Here, A s mn represents the symmetric anisotropy. Due to 
the space-group symmetries, all three types of magnetic 
couplings belonging to the eight planar Ti-Ti bonds may 
be obtained from those of a single bond, and so is the 
case for the four inter-planar bonds, see Tabic VII. 

The various magnetic couplings appearing in Eq. (18) 
are obtained by perturbation theory to leading order in 
V mn , namely, to second order in the hopping t mn and to 
first and second order in the spin-orbit coupling (scaled 
by A). In order to accomplish this calculation, we in- 
troduce the projection operator P mn onto the ground- 
state of H^ n , and the combined resolvent and projec- 
tion operator S mn onto the excited states. 24 In terms of 
these projection operators, the various terms appearing 
in Eq. (18) acquire the following structure. The Heisen- 
berg isotropic exchange, to leading order in the Ti-Ti 
hopping, is 

t o o _ pO rrtunn rrtun pO (^Q) 
mn J m J n mn mn mn mn mn' \ J 

The second term in Eq. (18) is the Dzyaloshinskii- 
Moriya antisymmetric anisotropic exchange interaction, 
which arises from second-order processes in the tun- 
nelling Hamiltonian, and first-order processes in the spin- 
orbit coupling, 

F) fS y S "1 = P° fftunc Iftunn rrso pO 
mn\ m nj mn mn mn mn mn mn mn 

-4- J-J so rrtun q rrtun pO 

' mn mn mn mn mn mn mn' 



where k — m,n and a,j3 = x, y, z. Since the Hamiltonian 
is invariant under time-reversal, there are no single-ion 



(20) 

In fact, there are additional terms in this order, 
in which there appear two Ti 2+ resolvents, e.g., 
fl^mn^mnCii^iifmiSi' These are smaller than 
the ones we keep, by an additional factor of ~ A c f/J7 c fr = 
0.059, where A c f = 0.208 eV is the gap between the 
ground state of the single-particle crystal field and the 
first excited state, see Table III. Following Ref. 25, we 
denote the vectors D mn , which refer to the microscopic 
single-bond couplings of the spins, as the Moriya vec- 
tors. The macroscopic antisymmetric anisotropic cou- 
plings between the sublattice magnetizations of the clas- 
sical ground state (discussed in the next section) are re- 
ferred to as the Dzyaloshinskii vectors. They are related 
to the Moriya vectors but are not necessarily the same. 

Finally, processes which are second-order in both the 
tunnelling and the spin-orbit interaction, yield 
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m mn n ' mn \ m n) mn mn ran mn mn mn mn mn mn 

_i_ pO ZT so C ZT so C rrtun c rrtun pO _■_ pO rrtun o rrtun o rj so o rj so pO 

* ran ran ran ran ran ran mn ran ran * ran ran ^ mn ran ^ mn mn 1 ^ mn mn mn 

i pO TJ so Q% f7 so pO rrtun o rrtun pO _■_ pO rrtun o rrtun pO f/ so £7 so pO 

' mn mn^mn^ mn* mn n ran ^ mn 1 * mn mn ' ^mn^mn ^ mn* 1 mn J -mn £1 mn'^mn- n mn mn' 

(21) 



These terms give rise to the symmetric anisotropies A^ n , 
as well as to corrections , of order A 2 , to the Moriya 
vectors. We have again omitted terms including two Ti 2+ 
resolvents. 

As was shown in Ref. 25, a systematic description of 
the magnetic anisotropies due to the spin-orbit interac- 
tion requires both the first and the second order pro- 
cesses in A. The technical reason being that the expec- 
tation value of the cross product in the second term of 
Eq. (18) is, in fact, also of order A, so that altogether the 
Dzyaloshinskii-Moriya interaction is at least second order 
in the spin-orbit coupling. As a result, although the anti- 
symmetric Dzyaloshinskii-Moriya interaction alone gives 
rise to spin-canting, when taken together with the sym- 
metric anisotropy, the system may, under specific condi- 
tions, still preserve rotational invariancc of the spins. 

The detailed calculation of the various terms appearing 
in Eqs. (19), (20), and (21) is lengthy, albeit straight- 
forward. More details are given in Appendix C. The 
values we obtain, using the parameters cited above, are 
listed in Table VIII. A comparison with spin-wave mea- 
surements is given at the end of the following section. 

IV. THE CLASSICAL GROUND STATE 

A. The magnetic order of the classical ground state 

The single-bond spin Hamiltonian, Eq. (18), is the ba- 
sis for the magnetic Hamiltonian, from which the mag- 
netic order of the classical ground state follows. To con- 
struct the latter, the entire Ti-lattice is decomposed into 
four sublattices. Namely, each magnetic unit cell includes 
four Ti ions, just as the crystallographic unit cell. The 
four sublattices are hence enumerated according to the 
numbers of the four Ti ions per unit cell shown in Fig. 1 
(sublattice i = 1 corresponds to Ti ion 1 and its Bravais 
translations, etc.). Assigning a fixed-magnitude magne- 
tization (per site) to each sublattice, Mj, one sums over 
all bonds which couple the four sublattices, to obtain the 
macroscopic magnetic Hamiltonian in the form 

H M = E [W M i + D S ( M * X M .) + M i T ij M j] > 

ij 

(22) 

where ij runs over the sublattice pairs 12, 13, 24, and 34 
of Fig. 1. This summation procedure gives rise to the 
macroscopic magnetic couplings: 1^ is the macroscopic 
isotropic coupling, are the Dzyaloshinskii vectors (to 



leading order in the spin-orbit coupling A), which are 
the macroscopic antisymmetric anisotropies, and are 
the macroscopic symmetric anisotropy tensors (of order 
A 2 ). The relations between those macroscopic couplings 
and the microscopic single-bond couplings are listed in 
Table IX. The inter-relations among the macroscopic 
couplings, which are dictated by the symmetries of the 
space-group, are contained in Table X. 

Although all four sublattice magnetizations are of 
equal magnitudes, their directions are all different. Those 
determine the magnetic structure of the classical ground 
state. In general, according to the space group Pbnm 
symmetries, there are four possibilities for the symme- 
try of sublattice magnetizations of the classical ground 
state, as listed in Table XL 2 In our case, the minimiza- 
tion of the magnetic Hamiltonian, Eq. (22), yields the 
magnetic structure shown in Fig. 3, drawn according 
to Table XII. This ground state has the following sym- 
metry: The x components of the magnetizations order 
antifcrromagnetically, in a G-type structure. The y com- 
ponents order antiferromagnetically as well, but in an 
A-type structure. Finally, the z components of the mag- 
netizations order ferromagnetically. Consequently, the 
classical ground state is characterized by two canting an- 
gles, ip and whose values are given in Table XIII, (see 
also Fig. 3). Due to the dominating Heisenberg coupling 
one observes that the magnetic structure of the classical 
ground state is predominantly G-type. The easy direc- 
tion along the x axis and the canting angles (both pro- 
portional to the spin-orbit coupling A) result from the 
anisotropic couplings of the model. Those break the ro- 
tational invariance of the magnetizations, and also cause 
the deviations from the pure G-type structure. 

Our magnetic structure is fully consistent with the ex- 
perimental one, as reported in Ref. 4. This experiment 
(in contrast to the one reported in Ref. 2) reveals that 
the G-type structure is indeed along the x direction, while 
the ferromagnetic moment is along the z direction. More- 
over, since the experiment of Ref. 4 is not sensitive to a 
small moment along the y axis, 26 our small A-type anti- 
ferromagnetic order along this direction does not contra- 
dict the data. We emphasize again that symmetry allows 
for such ordering, given the G-type order along x and 
the ferromagnetic order along z. Indeed, in the YTi03- 
system, which has the same space group as LaTiC-3 , such 
order has been detected, 27 but with different magnitudes 
of the canting angles, which cause the ferromagnetic or- 
der to dominate. 

One should note that by using naively the proce- 
dure outlined above to obtain the energy of the classical 
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magnetic ground-state, one obtains in the energy non- 
systematic contributions up to fourth order in the spin- 
orbit coupling A. To exemplify this point, we consider 

(H M ) = [A : ] -2(/ 12 + J 13 ) cos^cos 2 



the expectation value of i?M for the Ti-ions pairs r? = 12 
and 13, expressed in terms of the angles tp and $ and the 
superexchange couplings, 



[A 2 : ] +2(7 12 - J 13 ) sin 2 p cos 2 + 2(7 12 + J 13 ) sin 2 + 4(D 12 



Dy 



D® 3 V ) cos </? cos # sin $ 



+4£)? 2 2 cos tp sin p cos 2 - 2(Tff + r?f ) cos 2 p cos 2 
[A 3 : ] -4.0^ sin 93 cos t? sin 1? - 4rff cos <p sin <p cos 2 1? 

[A 4 : ] +2(r™ - T\l) sin 2 p cos 2 + 2(rf| + rff) sin 2 - 2ff| sin p cos sin 0. 



(23) 




FIG. 3. The magnetic order of the Ti ions in the classical 
ground state of the effective spin Hamiltonian of the lattice. 
The ions are enumerated according to the sublattice to which 
they belong. The x components of the spins order antiferro- 
magnetically in the G-type configuration, the y components 
order antiferromagnetically in the A-type one, and the z com- 
ponents order ferromagnetically. 

The leading orders of the terms are indicated in the 
square brackets. The non-systematic contributions of 
fourth order in A are due to the couplings r™, rf|, rff, 
rff, r| 2 (which arc all of order A 2 but are multiplied 
by sin sin 2 tp and sin $ sin y which are also of order 
A 2 ), and the A 2 correction of D^ 2 (which is multiplied by 
sin § sin tp) . Those contributions have been excluded from 



our calculation of the canting angles. On the other hand, 
we do include in the minimization of (i7 M ) terms up to 
the third order in A. This implies a systematic deriva- 
tion of the canting angles to first order in this coupling. 
(The classical ground-state energy has been found consis- 
tently, term by term, to second order in A. Although Tff 
and Tff have been calculated only up to order A 2 and 
consequently, we do not have the complete third-order 
term, this is of little importance when the canting angles 
are determined, since those couplings appear only with 
cost? and cosy, and therefore just cause an energy shift 
in (H M Y) Note also that although the Dzyaloshinskii 
and Moriya vectors first appear in linear order in A and 
the symmetric anisotropy coefficients in quadratic order, 
both kinds of anisotropies have to be considered as they 
cause terms which contribute in the same order of A to 
the classical ground-state energy. 

The values we obtain for the angles tp and 1? are listed 
in Table XIII. In particular, 1? = 0.80°, agrees, within 
the experimental error, with the values reported in Rcfs. 
2 and 11, 0.85° and 1.5°, respectively. (Note that 1? is 
defined here with respect to the oh plane.) This cant- 
ing angle causes the weak ferromagnetic moment, of or- 
der 0.014 ^b- This value agrees with the experimental 
one, within the uncertainty of the measurements which 
is caused by twinning of the crystal. 28 

Recently, an attempt has been made to analyze the 
relation between the anisotropy of the spin couplings 
and the paramagnetic susceptibility, which also has some 
anisotropy. 23 In this work the anisotropy of the spin cou- 
plings is taken into account by postulating phenomeno- 
logically an xyz model, which couples neighboring Ti 
spins and corresponds in our calculation to the coeffi- 
cients Af nn from Tabic VIII. A model susceptibility, 
which results from the xyz coupling via a molecular-field 
approximation (and from single-ion as well as from cova- 
lcncc effects), is calculated and then fitted with a couple 
of free parameters onto the measured susceptibility. As 
we have shown in this section, the antisymmetric and off- 
diagonal symmetric anisotropies — in the particular case 
of Eq. (23), components of the Dzyaloshinskii vectors 
and — can have at least the same conceptional im- 
portance for the magnetic properties of LaTiC>3 as the 
xyz anisotropics. Alternatively spoken, in general the 
xyz anisotropy is not the dominant anisotropy. This 
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basic argument is not restricted to the low-temperature 
case, which is accompanied by the magnetic order and 
which we investigate in the present paper, but refers also 
to the underlying spin couplings which influence high- 
temperature properties like the paramagnetic susceptibil- 
ity (whereas, e.g., the difference between the structural 
parameters of the low- and high-temperature case might 
correspond to more or less slight differences of the spin- 
coupling coefficients and of the orbital ground state) . In 
the way of an extension of Ref. 23, the question whether it 
is possible to include also the other than xyz anisotropics 
via free parameters in a model susceptibility and to com- 
pare this susceptibility to experiment, might be interest- 
ing. 



B. Comparison with spin-wave data 

The magnetic order in the classical ground state is the 
common starting point for a spin-wave calculation. In 
the case of the spin Hamiltonian pertaining to LaTiC>3, 
Eq. (18), one expects a rich spin- wave spectrum. This 
calculation is currently being undertaken, and will be 
presented elsewhere. 29 Nevertheless, our results above 
may be roughly compared with the existing spin-wave 
data. To this end, we ignore the antisymmetric and the 
symmetric anisotropics and hence assume a classical Neel 
state (for which the spin- wave spectrum is gapless). 

Inelastic neutron scattering has yielded the same value, 
J = 15.5 meV, for the single-bond Hcisenbcrg coupling 
for both the Ti-Ti bonds in the ab planes and those 
in-between the planes. 6 This value has been confirmed 
by the evaluation of Raman spectra. 30 Were we to av- 
erage our calculated values over the six bonds of each 
Ti ion, we would have obtained a value which is 32% 
higher. This rather modest discrepancy can be easily 
removed, by fine-tuning the model parameters. For ex- 
ample, by using the value A e ff = 6.6 eV (as estimated 
from an LDA+DMFT calculation based on the recent 
structural data 31 ), or by using a smaller value for the 
Slater-Koster parameters, V p d a = — 2.2eV (keeping the 
ratio between those parameters fixed, V p dir = 1.2 eV) in- 
stead of — 2.4eV, 5 or any other combined reduction of 
both of these parameters. Since a detailed comparison 
with the data requires the full spin-wave calculation, wc 
do not attempt here any fine-tuning of the model param- 
eters. For the purpose of the present paper, it suffices 
that the calculated Heisenberg couplings are consistent 
with the experimental value, within the uncertainties of 
our model parameters. 

Our calculation predicts somewhat different values for 
the in-plane Heisenberg coupling, J12, and the out-of- 
plane one, J13, yielding the ratio S = J13/J12 — 79%. 
Such an anisotropy may be detected by comparing with 
the spin-wave dispersion, e(q), at selected points, q = 
(0,0, 7r) and q = (tt,0,0), in the Brillouin zone of an 
effective cubic lattice of a unit lattice constant. When 



only this anisotropy is taken into account, then linear 
spin-wave theory gives 



e(q) = J12 y (2 + S) 2 - (cosq x + cosq y + S cosq z ) 2 . (24) 

With our calculated 5, wc find e(0, 0, n) /e(n, 0, 0) = 94 %, 
well within the experimental error bar of about 10 % for 
the spin-wave energies, from which the equality of the 
Heisenberg couplings on all bonds has been deduced. 32 
Hence, the (approximate) isotropy of the spin-wave spec- 
trum due to Ref. 6, which has been used as an argument 
to support the orbital-liquid state, 3 is consistent with our 
model. 

The calculation of Ref. 14 yielded a different value for 
the Heisenberg coupling ratio, 5 = 106%, i.e. a larger 
coupling along the c axis. This discrepancy can be traced 
back to our different crystal-field spectrum. In our case, 
the hopping amplitude between the crystal-field ground 
states on neighboring Ti ions is about 10% smaller for 
the bond 13 than for 12, see Table IV. This is a geometric 
effect which follows from the structural data. 4 

From a fit to the observed spin-wave gap, of or- 
der A = 3.3 meV, in conjunction with a spin model 
including solely antisymmetric anisotropics, a value of 
D = 1.1 meV has been deduced for the magnitude of 
the Moriya vectors. 6 We obtain higher magnitudes for 
the Moriya vectors. However, a full spin-wave expan- 
sion based on the Hamiltonian (18) indicates that the 
spin-wave gap is in fact dominated by the symmetric 
anisotropies rather than by the antisymmetric ones. 29 It 
is the canting of the ordered spins with respect to each 
other which is dominated by the Dzyaloshinskii vectors. 



V. SUMMARY 

We have presented a detailed analysis of the magnetic 
order pertaining to the LaTiC>3 system. The starting 
point of our calculation is the Ti-d orbital configuration 
which results from the static crystal field that includes 
the Jahn- Teller distortion, and which gives rise to orbital 
ordering as found in the experiment. 4 This orbital order- 
ing rules out the orbital-liquid picture 3 for LaTiC>3, which 
ignores the Jahn- Teller-like t^g splitting scheme and the 
resulting non-degenerate orbital ground state. 

Employing a perturbation expansion of this non- 
degenerate ground state in the effective hopping between 
neighboring Ti-ions, and in the on-site spin-orbit coup- 
ling, we have derived an effective spin Hamiltonian. It 
includes, in addition to the Heisenberg isotropic interac- 
tion between nearest-neighbor Ti-ions, the antisymmet- 
ric Dzyaloshinskii-Moriya coupling, and the symmetric 
anisotropic coupling. These three interactions conspire 
together to yield the magnetic order. We have found, by 
minimizing the magnetic energy of the classical ground 
state, that the magnetic order is primarily that of a G- 
typc antifcrromagnet, with the ordered moment along the 
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crystallographic a axis, accompanied by a weak ferromag- 
netic moment along the c axis. This configuration is in 
good agreement with the experimental findings. In addi- 
tion, we have found that there is a small A-type moment 
of the spin components along the b axis, which (although 
not yet detected in experiment) is allowed by the symme- 
try of the system. We find that the in-plane Hciscnberg 
coupling energy is about 27% higher than that pertain- 
ing to the coupling between a6-planes. By using these 
values in a spin-wave theory for the Heisenberg coup- 
lings, we show that both couplings are consistent with 
the isotropic spin-wave dispersion measured by inelastic 
neutron scattering. 6 

Our method seems to be particularly suitable to de- 
scribe the ferromagnetic Mott insulator YTi03 as well. 



Preliminary calculations (to be presented elsewhere) in- 
deed indicate ferromagnetic couplings in the a6-planes. 29 
Since the covalent pd hybridization in this system is as 
strong as in LaTiC>3, it will be of much interest to com- 
pare the classical magnetic Hamiltonians of the two sys- 
tems, and the ensuing spin-wave spectra. 
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APPENDIX A: THE EWALD SUMMATION 

The Madelung sum for the Coulomb potential in the 
point-charge model is given by 

1,71 

where r = (a;, y, z) is a point on the Ti ion No. 1, whose 
center is taken as the origin. In Eq. (Al), 1 are the 



Bravais translations, a n are the basis vectors of the unit 
cell, and q n are the corresponding point charges. The 
prime on the sum symbol indicates that the ion at the 
origin, 1 = a n = 0, is omitted. This sum converges very 
slowly. One therefore uses the Ewald summation, 33 (see 
also Ref. 34) where the sum (Al) is mapped onto two 
sums which converge much better, and which can be com- 
puted to high accuracy. Using the Ewald summation, the 
Madelung potential can be expressed in the form 



V(r) ^-eJ^T^ e"^ +igr V ^e"^ - e V — erfc (G\l + a„ - r|) + erf (Gr). (A2) 

g #o Vc9 n 1,7* 11 + a " r| r 



Here g are the vectors of the reciprocal Bravais lat- 
tice, whose basis vectors are (2n/a, 0, 0), (0, 2ir/b, 0), and 
(0, 0, 2-k/c), V c = abc is the volume of the unit cell, and 
G is a frequency cutoff. The value of V(r) is of course in- 
dependent of G. This cutoff is chosen such that the sum 
over the real-space lattice and the one over the reciprocal 
lattice can be stopped after about the same number of 
sites, when the required numerical precision is reached. 
In Eq. (A2), erf and erfc are the error functions 

2 f z 

erf (z) = 1 - crfc(z) = —= e~ e ' dt. (A3) 
V 71 " Jo 

The Ewald construction requires the neutrality condition 
$>n = 0, (A4) 

n 

which is fulfilled in our case. 

In order to find the spectrum and the eigenstates of the 
static crystal field, we have replaced the potential V"(r) 
by the pseudo-potential V ps {v), which is its Taylor expan- 
sion including the second and fourth orders in r. These 
are the Taylor orders which have non-trivial matrix ele- 
ments with respect to the d orbitals. 35 For instance, the 
second (fourth) Taylor order includes also terms which 



are proportional to xy (x 3 y,x 2 y 2 ). For the Taylor ex- 
pansion we use (r = |r|) 

-erf(Gr) = ^[l-i^ + i-L + ...j. (A5) 

The potential V^, s (r) is a harmonic function, invariant un- 
der inversion of the coordinates. The diagonalization of 
the matrix (7) V^> s (r)|7'), where 7 and 7' denote the or- 
thorhombic d orbitals, gives the results listed in Table III 
for the static crystal field. This calculation requires the 
second and fourth moments of the effective ionic radius, 
defined by 

J f 2 (r)r 2+n dr = (r n ), n = 2,4, (A6) 
where f(r) denotes the radial part of the ci-orbitals. 16 

APPENDIX B: THE HOPPING AMPLITUDES 
BETWEEN THE CRYSTAL-FIELD GROUND 
STATES 

As is mentioned in the text, the effective Ti— Ti hopping 
matrix elements between the crystal-field ground states 
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in the ab planes are of the same order of magnitude as 
those between planes, i.e., ~ I ^13 1 - This is a some- 
what surprising result in view of the fact that there is 
fcrro-orbital order in the planes and (canted) antiferro- 
orbital order between them. However, as we show here, 
in strictly cubic symmetry one has 1 1°2 1 = | * 13 1 - The de- 
viations from cubic symmetry cause the slight difference 
between these two hopping amplitudes. 

Let us hence consider the cubic case, and employ the 
coordinate system x',y',z in which the Ti sites 1 and 2 
are on the x' axis (see Fig. 1; these coordinates are ro- 
tated by 45° compared to the orthorhombic ones). The 
crystal-field ground states are now linear combinations of 
the three degenerate ti g orbitals 

\y'z), \x'z), \x'y'). (Bl) 

The hopping amplitudes are proportional to the over- 
lap of the two pertaining orbitals. Let us consider for 
simplicity the base-orbitals according to Eq. (3), i.e., 



|2) = 




\x' z) + I' 


|3) = 




\x'z) + 



x'y')), 

*V», (B2) 



where 1, 2 and 3 denote the relevant Ti-ions (see Fig. 1). 
Were we to find the direct hopping between the Ti sites, 
we would have obtained for the overlaps the result 

(1|2) = 1, (l|3> = -±, (B3) 

which reflects the ferro-orbital and nearly antifcrro- 
orbital order, respectively. However, the effective Ti— 
Ti hopping that we consider is mediated by the oxygens 
located between the Ti ions. Then, in strictly cubic sym- 
metry, for each pair of Ti ions, one of the three t2 g or- 
bitals cannot hybridize. This 'inactiveness' of one of the 
orbitals 3 is a direct consequence of the cubic symmetry, 
as is portrayed in Fig. 4, and is the source of peculiar 
hidden symmetries in the cubic Hamiltonian. 12 ' 13 In our 
example, the orbital \y' z) is inactive for the 12-bond, 
while for the bond 13 the inactive orbital is \x'y'), and 
consequently 

(1|2) I( - (x'z\ + (x'y'\)( \x'z) + \x'y')) = §, 
<l|3>-4( (y'z\-(x'z\)(-\y'z) + \x'z)) = -§, 

(B4) 

leading to = \t°^\. 



operators dt. (d kia ) which create (destroy) an electron 
in the crystal-field eigen-orbital i with spin component a, 
on the Ti ion located at site k. However, it is more con- 
venient to treat the two-electron states (which appear in 
the intermediate stages of the perturbation expansion) 
using the orthorhombic basis, Eq. (1). We denote the 
operators pertaining to this basis by cf (c fc7<T ), where 
7 enumerates the orthorhombic orbitals. The first part 
of this Appendix is devoted to the transformation of the 
Hamiltonian between the two schemes, and the diagonal- 
ization of the two-electron states. In the second part, we 
summarize the detailed expressions of the various terms 
resulting from the perturbation expansion. 

1. The Hamiltonian 

Denoting the matrix of the crystal-field Hamiltonian 
in the orthorhombic basis by V(k), we have 

HCt = EW*) c U<w (ci) 

£71 72 cr 

The matrices V(k) are real and symmetric. We next in- 
troduce the (unitary and real) matrix W(k) which diag- 
onalizes the crystal-field Hamiltonian, bringing it to the 
form 

ff cf = E^4«A^ (C2) 

kicr 

where Ei are the crystal-field eigenvalues, listed in Ta- 
ble III. These single-particle energies are shifted so that 
E = OcV, E x = 0.209 eV, etc. The relations between 
the operators d\ ia and c\ are hence given by 

i=E^(*) c L. C = E^( fc )4„ (ca) 

7 i 

such that 

W{k)V(k)W t (k) = E, (C4) 

with E = diagE'i. The diagonalizing matrix pertaining 
to site 1, W(l), is given in Table III. All other W{k) and 
V(k) follow from the symmetry properties of the unit 
cell, and are given by 



APPENDIX C: THE EXPLICIT CALCULATION 
OF THE EXCHANGE COUPLINGS 

Here we document the technical details of the pertur- 
bation calculation, that yields the effective spin Hamil- 
tonian. 

As explained in Sec. II, we consider a cluster of two 
nearest-neighbor Ti ions. The Hamiltonian of this clus- 
ter, given in Eqs. (5) and (6), is expressed in terms of the 
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(a) The hopping between the orbitals \y'z) 




(b) The hopping between the orbitals \x'y') 




(c) The hopping between the orbitals jx'z) 




FIG. 4. In cubic symmetry, there is one t2 g orbital per each 
Ti-O-Ti bond, which cannot participate in the hopping be- 
tween the Ti and the O ions because of the parity of the 0-2p 
orbitals. This is shown in panels (a)-(c), for each of the three 
t2 g orbitals, respectively, (a) The hopping between the or- 
bitals at the Ti sites 1 and 2 is not possible, since on 
the intermediate oxygen site there is no p orbital with a par- 
ity that would allow such hopping. For the same reason, an 
electron cannot hop from the p orbital |a:'), which is shown 
between the Ti sites 1 and 2, to any of the t2 g orbitals of 
the Ti sites. On the other hand, the hopping between the 
orbitals at the Ti sites 1 and 3, which is mediated by 
the orbital |y'^ on the intermediate oxygen site, is possible, 
(b) Analogously to the case (a), the hopping between the or- 
bitals |a;'y') at the Ti sites 1 and 3 is not permitted, because 
on the intermediate oxygen site there is no appropriate p or- 
bital. Likewise, an electron cannot hop from the p orbital 
which is shown between the Ti sites 1 and 3, to any of the 
t2 9 orbitals of the Ti sites. The (allowed) hopping between 
the orbitals \x'y lS } at the Ti sites 1 and 2 is mediated by the 

orbital on the intermediate oxygen site, (c) The hopping 
between the orbitals \x'z^ at the Ti sites 1 and 2 is medi- 
ated by the orbital on the intermediate oxygen site. The 
hopping between the orbitals |a;'z) at the Ti sites 1 and 3 is 
mediated by the orbital \x'\ on the intermediate oxygen site. 
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a. The Coulomb Hamiltonian 

The Coulomb Hamiltonian, in the orthorhombic basis, is 
given by 



2 / ^ 71727374 

T1T2T3 T4 



(C6) 



In order to specify its matrix elements in the 3d 2 sector 
(taken from Rcf. 19) we construct the triplet ("J^ 7 ) and 

the singlet ( 1 J ,7 7 7 ) two-particle states in the orthorhombic 
basis, 



-J- 77 (fc; aa 1 ) 

*F'(W) 



2 y kj<r k-y a 

-^P(k;aa'), 



2 ( C fc7<T C fe 7 '(T' C Ijct' C Ij'ct) 



k^a k^a 1 ' 



(C7) 



where we have used the direct matrix product a ~b®c, 



i.e., ci^j — bijCij- 



Altogether, there are 10 triplets and 15 singlets (in the 
a 1 = —a sector). Enumerating the triplet states in the 
following order, 



1) = |2) = |3) = K) - *r . I 5 ) = *r 



|6) = 1 7) = * 



24 
T ' 



M/ 21 , |9)=^ 4 , |10> = **\ 



(C8) 



the Coulomb Hamiltonian in the triplet sector becomes 



Urn = 



- 8B 
































A-5B 


3B 

























3B 


A-5B 3BV3 

























wVs 


A + B 
































A-hB 


35 


-3By/3 























3B 


A-5B 


-3BV3 























-3BV3 


-3W3 


A + B 
































A-8B 
































A — 5B 


-65 


























-6B 


A + 4B 



T- 



(C9) 



Here, A, B, and C are the Racah parameters, given by combinations of the Slater integrals F 2 and F4, 
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A = A -—F 4 , 



5 F 

4iT F4 ' 



35 



(CIO) 



The parameter Aq is determined such that upon diagonalizing simultaneously the Coulomb Hamiltonian and the 
crystal-field one, the lowest state has the energy U e g as explained in Sec. II B. Enumerating the singlets in the order 



|1> 
|6> 

111) 



*| 5 , |2) = *| 2 , |3) = *| 3 ,|4) 
*| 2 , |7) = *|\ |8) = *| 4 , |9) 
*| 3 , |12) -*| 3 , |13) 



*5 4 , |5) = * 



s > 



vl/ 13 



|14) 



*| 4 , 1 15) = * 24 



s > 



(Cll) 



the Coulomb Hamiltonian in the singlet sector becomes 



13 



Uc 



where 



Usi = 



A + AB + 3C AB + C 
AB + C 
3B + C 
3B + C 
C 




Usi 
U S 2 
U S3 



3B + C 



(C12) 



3B + C C 

A + AB + 3C B + C B + C AB + C 

B + C A + AB + 3C 3B + C 3B + C -By/6 

B + C 3B + C A + AB + 3C 3B + C By/6 

AB + C 3B + C 3B + C A + AB + 3C 

-BVE By/6 A + 2C 







U S 2 



A - 



AB 





2C 





A + B + 2C 
2BV3 





2BV3 
A + 2C 



U S3 



A 



B 
3B 

By/3 







2C 



A 



3B 
VB + 2C 
By/3 






-BV3 
By/3 
A + 3B + 2C 










A + B + 2C 
-3B 

-By/3 



A - 







-3B 
B + 2C 
-By/3 







-B^/3 
-By/3 
A + 3B + 2C 



(C13) 



b. The spin-orbit Hamiltonian 

Written in the orthorhombic basis, the spin-orbit Hamil- 
tonian is 



rrso _ ^ \ ^ ja a t „ 



(C14) 



where a takes the values x, y, and z, a a are the Pauli 
matrices, and the angular momentum matrices L", are 



L x = 



L y 









—i 








iy/3 








- 


-iy/3 





—i 


i 














i 











i 











0- 


iy/3 





—i 














iy/3 





—i 








i 











2i- 


















i 










-i 







-2i 











(C15) 



Transformed into the crystal-field eigenstates, the spin- 
orbit Hamiltonian takes the form 



(C16) 



with 



L?Ak) =Y J W il {k)L^ l W t 1 , v {k). (C17) 

77' 

Note that the relation L" , = — implies that 
L?Ak) = -LUk). 



c. The diagonalization of the two-electron states 

When there are two electrons on the same Ti ion (at site 
k), their state is described by d\ i<7 d\- a , . Using Eqs. (C3) 
and (C7), we re- write this state in terms of the singlet 
and triplet states, 



= W il {k)W jl ,{k){^J l {k;aa')8 11 , 

T 



+^l[n Y (k;aa') + ^'(k;aa')\(l- 5^)}, (C18) 



where we have omitted for brevity the summation nota- 
tions. Adopting the enumeration conventions Eqs. (C8) 
and (Cll), this state can be cast conveniently into the 
form 



14 



10 

4Xv = £<(Mj)*£(W) 
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(C21) 



+ J2<(k;ij)V s (k;™')- (C19) 

Here we have introduced the 10-dimensional vector wt, 
whose components are 



W i2 (k)W j5 (k), 

W t3 (k)W j5 (k), 
W a (k)W j4 (k), 
W l3 (k)W j2 (k), 
W H (k)W j5 (k), 
W a {k)W 3l {k), 
W i4 (k)W j2 (k), 
W n (k)W j2 (k), 
W i4 (k)W j3 (k), 

W a (k)W j5 {k)\ (C20) 



w T (k;ij) = ^ ^[W i5 (k)W j2 (k) 

W t5 (k)W j3 (k) 
W l4 {k)W n (k) 
W l2 (k)W j3 (k) 
W tb {k)W j4 {k) 
W tl {k)W 33 {k) 
W i2 (k)W j4 (k) 
W l2 {k)W 3l (k) 
W i3 {k)W j4 {k) 

W t5 (k)W 3l (k) 
and the 15-dimensional vector ws, 



ws(k;ij) 



W t5 (k)W 35 (k),W l2 (k)W 32 (k), 
W t3 (k)W 33 (k),W l4 (k)W 34 (k), 



for the entries fi = 1, 5, and 



ws(k;ij) = 



W l5 (k)W 32 (k) 

W l5 (k)W 3l (k) 
W i3 (k)W j4 (k) 
W l2 (k)W 3l {k) 
W lA {k)W 3l {k) 
W i5 (k)W j3 (k) 
W l2 (k)W 33 (k) 
W l3 (k)W 3l (k) 
W i5 (k)W 3i (k) 

W l2 (k)W j4 (k) 



+ W l2 (k)W 35 (k), 

+ W il (k)W j5 (k), 
+ W i4 (k)W j3 (k), 
+ W a (k)W j2 (k), 
+ W tl (k)W j4 (k), 
+ W l3 {k)W 3b {k), 
+ W i3 (k)W 32 (k), 
+ W a (k)W j3 (k), 
+ Wa(k)W jB (k), 

+ W i4 (k)W j2 (k) 



(C22) 



for the entries ^ = 6, 15. 

In order to obtain the two-states energies, we need 
to diagonalize simultaneously the Coulomb Hamiltonian 
and the crystal-field one. We have already written the 
Coulomb Hamiltonian matrix in terms of the triplets and 
the singlets. The next step is to express the crystal-field 
Hamiltonian in terms of those. Omitting for brevity the 
site index k, the crystal-field Hamiltonian matrix in the 
triplet sector is 



V T = 



V22+V55 


v 23 





-v 35 


v 24 





-v 45 


-Vis 





v 12 


v 23 


v 33 +v 55 





v 25 


v 34 


Vl 6 










v 13 








Vn + V 44 





-v 15 


-v 34 


-Vn 


^24 


-Via 


v 45 


-v 35 


v 25 





v 22 +v 33 





v 12 


v 34 


v 13 


-v 24 





v 24 


v 34 


-v 15 





v 44 +v 55 





v 25 





v 35 







Vl 6 


-v 34 


V12 





V11+V33 





-^23 


-V 14 


-v 35 


-v 45 





-v 12 


v 34 


v 25 





^22 + ^44 


v 14 


v 23 





-v 16 





v 24 


v 13 





-^23 


v 14 


V n +V 22 





v 25 







-v 13 


-v 24 


v 35 


-v 14 


v 23 





^33 + ^44 





v 12 


v 13 


v 45 





v 14 


-v 35 





v 25 





^11 + ^55 



(C23) 



We are now in position to find the resolvent operator in 
the triplet sector. Denoting by B the (unitary and real) 
matrix that diagonalizes the triplet part of H ci + H c , and 
by Ej, the corresponding eigenenergies, we have 

00') = -J2 *r(M, (fc; ™'\ (C24) 



One notes that since B k (fi,[i') = /i), the matrices 

Xj. satisfy 



(C26) 



Turning now to the singlet sector, we first find the 
crystal-field Hamiltonian matrix of the singlets, 



where 1/A£ is the resolvent operator 24 and 



(Ul = l 



B fc(M,Mi)£fe(Mi,M') 



E^ 



(C25) 



Vo = 



Vsi V S2 V S3 
Vs2 V S4 V S5 
v s 3 V S5 V s& 



= vl 



(C27) 



where 
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Vsi 





n 
u 


u 


n 
u 


n 
u 


V ^ V25 


o 


v 22 


o 


o 


o 


V ^ ^25 








2^33 




















2^44 

























V2V 25 


V2V 25 











V22 + V55 



^S2 



v^Vis 

V2V12 

a/2V 34 

x/2^34 

V^Fis V2F12 

Via V 15 



V S3 = 






^35 








V2V45 











V2V 23 








V2V24 





V2V 35 


V2V23 


V2T/ 13 








V2V 15 











\/2V4 5 


V2V24 


V2V 14 








V2V 13 











^23 


v 35 





^24 


v 45 



V n + V 55 V 25 
V S4 = V 33 + V 44 

V25 Vh + V22 



, V S5 



V45 V13 U35 K14 

v 13 V 45 V 24 Vu v 35 v 23 
v 24 v 13 v 23 v u 





' Vn + V 44 








v 34 




v 12 - 




1 


4)3 + V55 


^25 


Vis 


v 34 





v S(i = 




v 34 


V25 

^15 


V22 + V33 

v 12 


v 12 

Vn + V33 






v 34 








V 34 








v 44 + V55 


v 25 




. v 12 





v 34 





v 25 


V 22 + V44 _ 



(C28) 



Then we introduce the (15 x 15) matrix C that diago- 
nalizes H cf + H c of the singlets and the corresponding 
eigenenergies Eg. Analogously to Eq. (C25), it is conve- 
nient here to define as well 

X^,') = -t Ck ^"if^"^ (C29) 



;ji=i 



(C30) 



which satisfies 

x%( Jl ,v!) = x k s { P !,ti. 

Analogously to Eq. (C24) we have 

00') = -J2 ™')- (C3i) 



Collecting the results above, the intermediate two- 
particle states of the perturbation expansion are now 
given in the form 

— d f 

a kia a k 3 <j> 

10 

= Y, v%{k-,ij)X* fa, at/) 
nn'=i 

15 

+ E (C32) 



The final step involves transforming back and 
into the d operators. Consider, for example, *S?j, with 
fi' = 7172- Using Eqs. (C3) and (C7), we find (omitting 
the summation notations for brevity) 



*J l72 (fc;t7CT') 

W il7l (fc)W i27a (fc) - Wi 27l (fc)Wi l7a (fc) 



(C33) 



A similar calculation holds for the singlets. We therefore 
may write 

^4.4^ - ^(«';^2)4 1CT 4 2CT < ; (C34) 



with 



10 



Z k (ii';i 1 i 2 )= ^ w!^(k;ii')X^(p, fi')w^, (k;iii 2 ) 

15 

+ £ /*')<' (fc;»i»2). 
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(C35) 

We note that since the X's are symmetric, it follows that 
Zk(ii'; = Z k (iii 2 ;ii'). (C36) 
Also, since Wg(k;ii') = u>g(k;i'i) and w^,{k;ii') = 



-Wj,(k; i'i), one has 



Z k (ii';iii 2 ) = Z k (i'i:i 2 h). 



(C37) 



2. Perturbation expansion 

Our formal expressions of the various terms in the per- 
turbation expansion, Eqs. (19), (20), and (21) above, 
involve the projection operators P^ n and S mn = (1 — 
P^ n )/A£. Here we give their explicit expressions in 
terms of the quantities derived in the first part of this 
Appendix. 

The projector onto the unperturbed ground-state 
space is 



Pmn — d m0a d n0a , \ o)(o\d n0a td 



mOa ! 



(C38) 



where |0) denotes the vacuum state. Similarly, the pro- 
jector onto the Ti 3+ sector of the excited states is 

PL = E d UAi 2 *> \0)(0\d ni2a ,d miia . (C39) 



era-' 
«1« 2 (^00) 



It follows that the resolvent operator applied to is 
given by 



act' 
«lt 2 (^00) 



(C40) 



In a similar way, the projector onto the Ti 2+ sector of 
the excited states is 

4=jE [E^(W)|0)(0|^(fc;aa') 

15 

+ E^(fc;^')|0)(0|^ t (fc;aa')], 

(C41) 

which gives, upon applying the resolvent operator, 

■^pL = -\H [E x T(M,M')^'(fc;^')|o) 



kaa' fifi' 



x{0\^(k;aa') 



x(0|^ t (fc;crcr') . 

(C42) 

Hence, the combined resolvent and projection operator 
onto the excited states is 



Q = - (p 1 4. p 2 "i 

u mn ^£ \ ran ~ ran)' 



(C43) 



Collecting these results, and expressing the products 
of d in terms of the ground-state spin operators [see 
Eqs. (15)], one obtains the magnetic exchange couplings. 
These are listed below. 



a. The Heisenberg couplings. These are given by 

J mn = - E 0i 2 )#£ + t™Z n (hO; Oia)*^] • (C44) 

b. TTie Moriya vectors. These are given to first order in the spin-orbit coupling, 

£>«„ = 2»A E ^{ifoM[C°^ m (n i ;0 J2 )e i +tj» Z n (i 1 0;0i 2 )Cn] 

[CSXM; oi 2 )c, + « 2 )ej }• (C45) 

c. T/ie symmetric anisotropics, and the A 2 correction of the Moriya vectors. These terms are of second order in A, 
and have a more complicated structure. In order to present them in a concised fashion, we write the left-hand-side of 
Eq. (21) in the form 

S m Al nn S n + 'D' mn (S m X S„) = J^ n S m S„ + E E\I7., ^mn(*)*')0^mn(0! (C46) 

where J' mn and D^ n are the A 2 corrections of the Heisenberg couplings and the Moriya vectors, respectively. In Eq. 
(C46) I enumerates the 4 spin invariants, such that 
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jafi (-1 \ _ pa a/3 r a/3 rt\ _ ™ o/3 , c/3 o« _ c o q 
Jmn^i — °m D n' 1 mn\ A ) ~ °m D n + D rn D n V ™ «' 



(C47) 



where e Q) 3 7 is the totally antisymmetric tensor. It remains to list the coefficients appearing in Eq. (C46). These are 
given by 



C^(*.i',2) 



mnj 



C^(i, i', 1) - 2A 2 £ Lf (m)i& (n) [&SX(«i; Oi 2 )t^ + Z^i'h; Oi 2 )C„ 

-C!2 ra (n0; i 2 i)<*! - t^Z n (hO; i 2 i')C, 
x2 E{^"oH4'N [*ji^m(i'ii;i2i)t^ + CXM; 0i 2 )t* 

+ ^«>(n)*&(") [4Um(*iO; Oi 2 )<M + t^^^OC] 
-L«(m)L Oi ,(n) [t^Z^O; ti 2 )t^ + C° Z n (iiO; i'i 2 )C„] }, 

+Lf (m)Lf, i (m) [C° «i0; *'* 2 )C 2 i + CXM; Ot 2 )tj£] }, 

A ' E {^M^w (») [&iXM; « 2 )t^ + Cm^n(^i; ^o)C„] 

+L? (n)Z&(n) 0t 2 )t^ + tj£Z n (iiO; t'ia)*^] }. 



C£&(m',3) 



C«i(i,i'A) 



(C48) 



TABLE I. The parametrization of the unit cell (space 
group Pbnm), modulo the lattice constants a,b,c. 



La 


(s;re,j/re, 1/4), (1/2 - xre, 1/2 + yuE , 1/4), 

(-IRE, -2/RE, 3/4), (1/2 + IRE, 1/2 - |/RE , 3/4) 


Ti 


(0, 1/2, 0), (1/2, 0, 0), (0, 1/2, 1/2), (1/2, 0, 1/2) 


Ol 


(a;oi,2/oi, 1/4), (1/2 - xoi, 1/2 + j/oi, 1/4), 
(-xou -Voi, 3/4), (1/2 + a;oi, 1/2 - j/oi, 3/4) 


02 


(xo2, yc>2, 202), (X02, j/02, 1/2 — 202), 

{-X02, -J/02, -202), (-£02, -J/02, 1/2 + 202), 

(1/2 - x 2, 1/2 + yo2, 202), 

(1/2 - x 02 , 1/2 + yo2 , 1/2 - 202), 

(1/2 + x 2, 1/2 - 2/02, -202), 

(1/2 + Z02, 1/2 - 3/02, 1/2 + 202) 



TABLE II. The symmetries of the space group Pbnm. 



Inversion centers 


Ti sites and centers of ab-planar Ti 
plaquettes, i. e. (0, 0, 0), (1/2, 1/2, 0), etc. 


Mirror planes 


2 = ±1/4 


Glide planes 


x = ±1/4, translation by (0, 1/2,0) 


Screw axes 


Through the inversion centers, along the 
2 axis, rotation by 180° 



TABLE III. The static crystal field for Ti 3+ (site 1): Spec- 
trum and eigenstates in the orthorhombic basis for the d or- 
bitals, see Eq. (1) and the following comment there. [The 
eigenenergies Ei and the matrix W(l) used in conjunction 
with Eq. (C4) are defined by the spectrum and the coordi- 
nates of the eigenstates, respectively, as given in this Table, 
where the first row of W(l) is the coordinate vector of the 
ground state, etc.] 



-0.468 eV 


(-0.035, 0.016, 0.770, -0.035, 0.636) 


-0.259 eV 


(-0.052, -0.397, 0.088, 0.911, -0.049) 


-0.239 eV 


(-0.407, 0.035, -0.587, 0.086, 0.693) 


0.452 eV 


( 0.853, 0.315, -0.197, 0.221, 0.290) 


0.515 eV 


(-0.319, 0.861, 0.123, 0.336, -0.169) 


Basis 


xy), 2z 2 ), yz), xz), x 2 ~y 2 ) 
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TABLE IV. The effective Ti-Ti hopping matrices for the 
d eigen-orbitals of the crystal field from Table III; values are 
given in eV. The rows and the columns are ordered beginning 
with the ground state of the crystal field (index 0), continuing 
with the first excited state (index 1), etc. t\3 is symmetric 
due to the mirror plane at z = 1/4. As the matrices are given 
in terms of the crystal-field eigen-basis, the dependence of 
the hopping matrices on the bonds is particularly simple. [In 
contrast, had we used the orthorhombic basis (1), additional 
minus signs would have appeared in several entries.] 



Planar 



— *25 — 


^65 — *34 


— f* - 

— z 38 — 


*47 — ^87 




-0.198 


-0.155 


-0.052 


-0.022 


0.016 


0.109 


0.133 


0.022 


-0.089 


0.135 


-0.114 


0.167 


-0.188 


-0.098 


0.193 


-0.021 


0.088 


-0.235 


0.579 - 


-0.710 


0.010 


-0.019 


-0.003 


0.089 - 


-0.121 



Inter-planar 



— *39 — 


^410 








0.178 


0.047 


-0.143 


0.010 


0.020 


0.047 


0.244 


0.072 


0.135 


0.224 


-0.143 


0.072 


0.146 


-0.008 


-0.057 


0.010 


0.135 


-0.008 


-0.112 


-0.312 


0.020 


0.224 


-0.057 


-0.312 


-0.812 



TABLE V. Model parameters of the calculation. 



Momenta of the effective ionic radius for Ti 3+ 
(r 2 ) = 0.530 A 2 , (r 4 ) = 0.554 A 4 
Slater integrals for Ti 2+ 
F 2 = 8F 4 /5 = 8.3 eV 
Spin-orbit parameter 



A = 18 meV 



Slater-Koster parameters 



V pda = -2.4 eV, V pd * = 1.3 eV 
Effective charge-transfer energies (Ti-Ti, Ti-O) 
U eB = 3.5 eV, A eff = 5.5 eV 



TABLE VI. The combined static and covalent crystal fields 
for Ti 3+ (site 1): Spect rum and eigenstates in the orthorhom- 
bic basis for the d orbitals. The covalent contribution is cal- 
culated for a TiC>6 cluster. The full eigenstates are linear 
combinations of Ti 3+ states and Ti 2+ states (accompanied by 
a p hole on one of the oxygen sites) . Here, only the Ti 3+ parts 
of the five lowest eigenstates are shown, corresponding to the 
states jd 1 ) of Eq. (13). 



-0.663 eV 


(-0.029, 0.020, 0.778, -0.023, 0.627) 


-0.441 eV 


( 0.084, -0.383, 0.200, 0.875, -0.200) 


-0.430 eV 


(-0.393, -0.089, -0.550, 0.279, 0.676) 


0.737 eV 


( 0.856, 0.322, -0.211, 0.173, 0.297) 


0.797 eV 


(-0.323, 0.861, 0.092, 0.354, -0.144) 


Basis 


xy), 2z 2 ), \yz), \xz), :r 2 -y 2 ) 



TABLE VII. Symmetries of the effective spin Hamil- 
tonian due to the space group. The relations among 
the anisotropic couplings are abbreviated as follows. 
(+,+,+)i2 = (-,+,+)ie means D 12 = (-Df 6 , D y 16 , Df 6 ), 
etc. Due to the mirror plane z — 1/4, the inter-plane Moriya 
vectors have vanishing z components and the inter-plane sym- 
metric anisotropies have vanishing yz and xz entries. The 
transformation of the symmetric anisotropies is characterized 
by the off-diagonal coefficients (A^ n , A^ n ,A^ n ) whereas the 
diagonal coefficients are invariant in and between the planes, 
respectively. 



Heisenberg couplings 


J12 = J16 = 


= J25 = Jes — ^34 = J38 = 


J47 = J&7, 




J 13 = J24 = J39 = J410 




Moriya vectors 


(+,+,+)l2 = 


(-,+,+)l6 = (+,-,~)25 


= (-, -, -)6B 


= +)34 = 


(+,-,+)38 = (-,+,-)47 


= ( + ,+,-) 8 7, 


( + ,+,0)l 3 = 


(+,-,0)24 = (-,-,0)39 = 


= (-,+,0)410 


Symmetric anisotropies 


(+,+,+)l2 = 


( + , -, -)l6 = (+, -, -)25 


= ( + ,+,+)65 


= (-,-,+) 3 4 = 


( — , +, — )38 = ( — , +, —)47 


= ( — , — , +)87, 


(0,0,+) 13 


= (0,0, -) 24 = (0,0,+)39 = 


= (0,0,-) 4 io 
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TABLE VIII. The calculated single-bond spin couplings 
(in meV). The Moriya vectors are given including the cor- 
rections Djnn, which are of order A 2 . The symmetric 
anisotropies are given as A^ n = (A^ n ,A]™ n ,A™ n ) and 
A^ d n = (CnX»i^2i) for the diagonal and off-diagonal 
entries, respectively. 

Heisenberg couplings 
J12 = 24.616, J13 = 19.416 
Moriya vectors 
D12 = (3.254, -1.273, -1.286), D i3 = (-2.886, 0.543, 0) 
Symmetric anisotropies 
A? 2 = (0.188, 0.066, 0.037), A d 3 = (-0.039, -0.218, -0.190), 
A°2 = (-0.035, -0.111, -0.088), A° 3 = (0, 0, -0.074) 



TABLE IX. The macroscopic couplings of the sublatticc 
magnetizations in terms of the microscopic single-bond spin 
couplings. For instance, we have ii 2 = J12 but /13 = Ji 3 /2. 
This is because the coordination number of a Ti ion is 4 in 
the planes and 2 between the planes. 

Isotropic couplings 
I12 = J12, Iva = \Jl3 
Dzyaloshinskii vectors 

D?2 = (0,-P|' 2 ,-Pl 2 ),Dl3 = |Dl3 

Macroscopic symmetric anisotropies 

r? 2 = a? 2 , r? 2 d = (as, 0,0), r 13 = \a 1z 



TABLE X. Symmetries of the magnetic Hamiltonian due 
to the space group. The relations for the anisotropic couplings 
are abbreviated as in Table VII. Due to the glide planes, the 
Dzyaloshinskii vectors of the planar bonds have vanishing x 
components, and the respective symmetric anisotropies have 
vanishing xz and xy entries. Because of the mirror planes, the 
Dzyaloshinskii vectors of the inter-planar bonds have vanish- 
ing 2 components and the respective symmetric anisotropies 
have vanishing yz and xz entries. 

Isotropic couplings 
I12 = I34, Irs = I24 
Dzyaloshinskii vectors 

(0, +, +)l2 = (0, -, +)34, ( + , +, 0)l3 = ( + , -, 0)24 

Macroscopic symmetric anisotropies 



TABLE XL All types of magnetic order which are allowed 
by the space group Pbnm. There are four possibilities, de- 
noted by x s , x 3 ", z s , and z a . They are allowed because the 
ordered state can be symmetric or antisymmetric according 
to the glide plane x = 1/4 and the mirror plane z = 1/4, 
respectively. The order in LaTiOs is of the first type. Here, 
G x denotes G-type antiferromagnetic moment along x, A y de- 
notes A-type antiferromagnetic moment along y, and F z de- 
notes ferromagnetic moment along z. The other possibilities 
involve also C-type ordering, e. g. C z for the z components of 
the magnetic moments. The magnetizations Mt of the sub- 
lattices are given in terms of Mi. (+, +, +)i = ( — , +, +) 2 
means M 1 = (-Mf , Mf , M 2 Z ), etc. 



1. X s , 2 a 


2. x% 2 a 


3. x s , z s 


4. x a , z s 


(+, +, +)l 


(+, +, +)l 


(+,+, +)i 


(+, +, +)i 


= (-, +, +)2 


=(+. -. -h 


=(-, +, +)a 


=(+, -, -h 


= (". +)s 


=(- -, +)a 


=(+, +, ~)s 


=(+, +, -)s 


= (+, +)4 


=(-. +. ~)4 


= (". +. -)4 


= (+, -, +)4 


Gx Ay Fz 


A X Gy C Z 


Cx Fy A z 


F x Cy G Z 



TABLE XII. The structure of the magnetic order (the first 
possibility of Table XI) , characterized by the sublattice mag- 
netizations Ms, in the classical ground state, in terms of the 
canting angles ip and 1?. Each of these angles is proportional 
to A, the spin-orbit parameter. 

x components: G-type 
-Mf = Mf = Mf = -Mf = cos ip cos ■& 

y components: A-type 
-Mf = -Mf =M% = Ml = sin ip cos ■& 
z components: ferromagnetic 
M{ = Mf = Mf = Mf = sint? 



TABLE XIII. The values of the macroscopic magnetic cou- 
plings in meV, the resulting canting angles of the magnetiza- 
tions in the classical ground state, and the resulting abso- 
lute values of the ordered moments (normalized to 1 /jb, i- e., 
without the reduction due to quantum fluctuations and the 
spin-orbit coupling). The three coefficients of the macroscopic 
symmetric anisotropies are taken into account (see text). 

Isotropic couplings 
7i2 = 24.616, /is = 9.708 



Dzyaloshinskii vectors 
D? 2 = (0, -1.273, -1.286), D? 3 = (-1.589, 0.271, 0) 
Macroscopic symmetric anisotropies 

rff = 0.188, rff = -0.020, r^f = -0.037 

Canting angles 
p = 1.42°, ■& = 0.80° 
Ordered moments 
M = (0.9996, 0.0248, 0.0140) mb 



(+, 0, 0)12 = (-, 0, 0) 34 , (0, 0, +)i 3 = (0, 0, -) 24 
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